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^ ■ Abstract 

, We study the scattering of the quantized electromagnetic field from a linear, 

' dispersive dielectric using the scattering formalism for quantum fields. The 



00 

' medium is modeled as a collection of harmonic oscillators with a number of 

o 

. sion, which is widely used to describe experimental data for real dispersive 

, media. The integral equation for the interpolating field in terms of the in field 

is solved and the solution used to find the out field. The relation between the 
in and out creation and annihilation operators is found which allows one to 
P5 ■ calculate the S-matrix for this system. In this model, we find that there are 

^ ■ absorption bands, but the input-output relations are completely unitary. No 

Q^. additional quantum noise terms are required. 

> ' 

X ■ I. INTRODUCTION 



A fundamental problem in quantum optics is how the properties of light change as it 
propagates through a medium. If the medium is nonlinear, new frequencies can be produced 
and the quantum noise properties of the field can be altered. This leads to such interesting 
phenomena as solitons, squeezing, and quantum phase diffusion, all of which have been ob- 
served If the medium is linear, the situation is not as dramatic, but linear media serve 
as a first step in the description of their nonlinear brethren, and they present problems in 
their own right, such as the inclusion of dispersion. It is generally thought that an accurate, 
first-principles treatment of dispersion necessitates the inclusion of absorption, and conse- 
quently additional noise or reservoir operators. Here we analyse a quantum field theoretic 
model that demonstrates dispersion-induced absorption, but without any additional noise 
operators appearing in the scattering relations. 

Dielectric media can be described in a number of different ways. They can be charac- 
terized by their susceptibilities, a procedure which, when fields are included leads to the 
macroscopic Maxwell equations. Consequently, we shall call this the macroscopic approach. 
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It leads to difficulties when one wants to include dispersion in a Lagrangian or Hamilto- 
nian formulation, because dispersion is a consequence of the fact that the response of the 
medium is not instantaneous, but depends on the values of the field over a range of times 
0. For certain kinds of fields, in particular narrow-band ones, these problems can be over- 
come by using an approximate Lagrangian which is local in time [0]. Another approach is 
to construct a microscopic model for the medium and to include the degrees of freedom of 
the medium in the theory This, for obvious reasons, we shall call the microscopic 

approach. It has the advantage that the inclusion of dispersion is not a problem, but the 
disadvantage is that for each new medium, a new model must be constructed. There are 
also intermediate approaches which use frequency-dependent susceptibilities, but also add 
quantum noise operators to the equations of motion for the fields 0, |1TT| . 

Most of the work on quantized fields in media has concentrated on what happens inside 
the medium. However, there has been a steady stream of research which has considered 
fields entering and leaving a medium as well. This is essential if one wants to describe real 
experiments, in which the fields are generated outside the medium, then pass through it, 
and are finally measured in free space. 

Perhaps the first to investigate this question were Lang, et. al., who examined the con- 
nection between the field inside and outside a laser cavity while studying why the laser 



linewidth is so narrow [12], |13]. Their model consisted of a cavity bounded on one end by 



a perfectly reflecting wall and on the other by a thin dielectric slab, and the cavity itself 
is filled with an active medium. This cavity is embedded in a larger cavity which repre- 
sents the universe. They showed how the modes of the universe are related to the cavity 
quasi-modes, which makes it possible to find the output field in terms of the cavity field. In 
their analysis the field was classical, but soon thereafter the quantum version of their model 
was constructed and used to investigate the relation between the field inside and outside 
the cavity for a laser in the linear regime by Ujihara fl^. This approach was later used by 
Gea-Banacloche, et. al. to study the relationship between the squeezing generated inside a 
cavity to that outside the cavity |T^. 

This last problem had ffist been considered by Yurke who based his approach on an 
earlier paper by Denker and himself They developed a quantum theory of electronic 
networks in which the network itself is located at a; = and transmission lines extending 
from there to x = +oo bring input signals to the network and carry output signals away. 
Fields propagating toward x = are input fields and those propagating away are output 
fields, and the object is to find the output fields in terms of the input ones for a given 



network. A related input-output theory was developed by Collett and Gardiner [|T^ and 
was put on a ffimer footing by Carmichael |T^. This theory considers a cavity containing a 
medium, active or passive, linear or nonlinear, which is coupled to a reservoir. The reservoir 
operators serve as the input and output fields. The dynamics of the system inside the 
cavity is described by a master equation, and its solution is used to find the time- dependent 
reservoir operators and, thereby, the output field. 

More recently a number of groups have examined the scattering of the quantized electro- 
magnetic field from inhomogeneous linear dielectrics. Glauber and Lewenstein considered 
the case of a non-dispersive, lossless dielectric which is described by a real position-dependent 
susceptibility |jl9[. The scattering from dispersive media was studied by Knoll and Leonhardt 
who considered a medium consisting of damped harmonic oscillators M. Rather than use 
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a formal scattering approach they used time-dependent Green's functions to solve the field 
equations. A different treatment of dispersive media was given by Matloob, et. al. They 
considered an arbitrary complex, frequency- dependent dielectric function and quantized the 
theory at the level of the equations of motion rather than starting with a Lagrangian. Work- 
ing if frequency space they found the fields emerging from a dielectric slab in terms of those 



entering it. A similar analysis was carried out by Gruner and Welsch |11|. A final approach 



is based on polaritons in finite media [^- [^. For an infinite dielectric interacting with 
the electromagnetic field the eigenstates of the Hamiltonian are mixed matter-field modes 
known as polaritons. If the medium is finite, the polaritons acquire a finite lifetime. By 
looking at the the electromagnetic parts of the polariton modes, scattering of the field from 
the medium can be described. 

In this paper we shall apply the quantum scattering formalism for fields to describe an 
electromagnetic wave scattering from a finite medium, which behaves as a mirror or beam- 
splitter. The medium is treated microscopically, and it is dispersive. The final result is an 
explicit expression for the out operators in terms of those of the input field. The calculation 
starts from first principles, and, consequently, shows how some of the relations between 
in and out operators, which are often used in quantum optics, follow from an underlying 
scattering theory. An important feature of the model used here is that it includes a dielectric 
medium with multiple bare resonances, leading to a number of discrete absorption bands. 
This is typical of real dielectric materials, and leads to a dielectric constant that can be 
modeled with the widely used Sellmeir [jlO[ expansion in frequency. 

As we noted in our discussion of previous results, there are three treatments of scattering 
from a dispersive, linear dielectric. Two are to some extent phenomenological in that they 



do not start from a Hamiltonian describing the field- medium system [^], [jTT|. The third 
approach treats time-dependent fields rather than finding the asymptotic in and out fields 
which are the basic objects in a scattering treatment [Q. We believe that this leaves room 
for a more fundamental approach which can place the theory on a firmer foundation. The 
results we find from field theoretic scattering theory are similar to those found in the ap- 
proach pioneered by Yurke, and have the virtue that they are simple and intuitively clear. 
By employing a fundamental approach, we have the advantage that the meaning of all of 
the operators which we employ is well-defined, which is not always the case in the more 
phenomenological treatments. The results presented here can be viewed as a justification of 
earlier phenomenological theories. 



II. MODEL 

We shall consider a one- dimensional model of the electromagnetic field and the medium 
which was developed in reference p|. This model can be used to describe the normal 
incidence of an electromagnetic wave on a medium, where the wave travels in the x direction 
and is polarized in the z direction. 

The field can be represented by means of the dual potential, A(x, t), which is appropriate 
if there are no free charges. In the case of a 2;-polarized normally incident plane wave, A(x, t) 
is the y component of the dual potential. The fields are given by 
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The medium consists of dipoles which are harmonic oscillators with masses rrii, and bare 
frequencies Q^, where u = 1, . . . ,N. The un-renormalized oscillator frequencies can be 
chosen to correspond to transition frequencies of atoms or molecules making up an actual 
material. Each oscillator is described by a field, r,^(x), which gives the displacement of the 
oscillator at position x and with frequency Q^. It is convenient to represent the oscillators 
in terms of the polarization fields 



Pu{x) = qvPu{x)r^{x) 



(2.2) 



where pu{x) is the density of oscillators with frequency Q^, and the dipole corresponding to 
oscillators of type v consists of charges q,^. We shall work in the multi-polar gauge so that the 
coupling between the electromagnetic field and the medium is proportional to Y.uPy{x)D{x). 
The medium self-interaction terms proportional to the square of the total polarization are 
incorporated into the frequencies, Vty. 

For a volume such as a waveguide of cross-sectional area A, the Lagrangian density for 
the medium-field system is given by 



A 



C = — \-K\x)-{d^k{x)f + Y: 
+2p,{x)d^k{x)]} , 



1 



where 



guix) 



(2.3) 



(2.4) 



A. Refractive index 



^From the Lagrangian density we find the equations of motion for the fields 



dtPu + ^Ipu = 9ud^A . 



(2.5) 



For a medium of constant density, (i. e. g^ix) is independent of x), we can solve the above 
equations by assuming that both p„ and A are proportional to e*^'^^"'^*) . The values of u are 



the frequencies of the modes of the system and are given by the solutions of equation |p.5 



{kc) 



1-E 



9v 



Defining the index of refraction, n{uj), to be kc/u, we find 

9u 



(2.6) 



n{Lu) 



.-1/2 



E 



(2.7) 



This is very similar [§] to the classical Sellmeir expansion for the refractive index. Note 
that this expansion is not identical to the Sellmeir expansion, but can be converted into the 
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commonly used Scllmeir form through a renormahzation of the bare resonant frequencies 
of the oscillators. The characteristic property of this type of equation is that it possesses 
solutions for the refractive index that are either purely real (transmission bands) or purely 
imaginary (absorption bands). At the bare resonance frequency, the refractive index is zero. 
Near a resonance, where a; — > fi;^, the refractive index is real for a; > Qj,, and imaginary for 
u! < At a finite detuning below a resonance, the refractive index goes to infinity just 
below the start of the corresponding absorption band. 



B. Lagrangian quantization 

^Prom the Lagrangian density we can find the canonical momenta corresponding to A 
and pu, which we shall denote by U and tTj/, respectively. These are given by 

U{x) = i^oHx) Mx) = . (2.8) 

eoQAx) 

The theory is quantized by imposing the commutation relations 

[A{x,t),fl{x',t)] ^ih6{x-x')/A (2.9) 

and 

[p,{x, t),7r,ix', t)] = ihS,yS{x - x')/A . (2.10) 

The canonical momenta and the Lagrangian density can now be used to find the Hamiltonian 
density for the quantized theory 

nix) = A : S^fl^x) + (dM^))' + Y:[elg.ix)nlix) 

+^Mx) - 2p,{x)dMx)\] ■■ ■ (2.11) 
9 Ax) J 

This can be put in a difi^erent form if we define annihilation and creation operators, iu{x) 
and Q{x), for the oscillators, where 



so that 

[Ux),il{x')]=6,y5{x-x')/A. (2.13) 
We finally have for the Hamiltonian density 

(2.14) 



A 
2^ 


: i^fl'ix 


-'\ 




1 ^ 
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III. SCATTERING THEORY 



In order to determine what happens when an electromagnetic wave scatters off of the 
medium we shall apply the standard formulation of scattering for quantum fields p3[. This 
is done in the Heisenberg picture so that it is the field operators which are time dependent. 
Because we shall consider a medium which is bounded in the x direction, the interaction 
is bounded in time. This can be seen either by considering the incoming waves to be wave 
packets, so that the interaction takes place only while the packet is inside the medium, or by 
using plane waves and turning the interaction on and off adiabatically. In either approach, 
the fields will go to free fields both as t ^ — oo and as t oo. The free fields as t — cxd 
are the in fields, and those as t ^ oo are the out fields. The time dependent field operators 
which carry the full time dependence of the Hamiltonian, including the interaction, are also 
known as the interpolating fields, because they interpolate between the in and the out fields. 
Our goal is to use the interpolating fields to find an expression for the out fields in terms 
of the in fields. This will give us a complete description of the scattering process. We note 
here that a related Heisenberg-picture approach to quantum scattering theory relevant to 
quantum optics measurements was recently developed by Dalton et al. They developed 
a similar basic formalism, but did not consider specific examples. 



To find the relationship between the in and out-fields, we need to express the equations 
of motion of the interpolating fields as integral equations. From the Hamiltonian for our 
model we find 



A. In and out fields 




(3.1) 



In order to express these as integral equations we define the Green's functions A^'^'^^\x,t), 
A('^'^'')(a;,t), T^;^^\x,t), and r(f'^'')(a;, t). They satisfy the equations 



{d^ - c'dl)A<^^^^'\x,t) = 5ix)6{t) 
{dt + tn,)T^J''\x,t)=6{x)6{t) 
{d^ - c^dl)A'-'"'^\x,t) = 6{x)6{t) 



{dt + zn,)T(j''\x,t)=6{x)6{t) , (3.2) 



and the boundary conditions 



A^'-^'\x,t) = T'i^''\x,t) 
AH^)(x,t) = r[^'^^\x,t) 



fort < 



fort > . 



(3.3) 



The retarded Green's functions can be expressed as 
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^ ' ^ (27r)2 

r(-)(x,t) = ^ 



dkduo 



— ie ~ uj) 



(3.4) 



where e — > 0^, and the advanced Green's functions are given by almost identical expressions, 
the only difference being that e is replaced by — e. The integral equations corresponding to 
the differential equations, Eqs. (|3J]) , are 



A(x, t) = A,„(x, t)-c^ I dx' I dt'A("^*)(x -x',t- t') 



Iheoguix') r 



^^{x, t) = t) + ^Jdx' J dt'T^'''\x -x\t- t') 



2n„ 



eoh 
d.^>A{x',t' 



(3.5) 



The corresponding expression involving the out-fields is: 

A(x, t) = Aoutix, t)-c'' f dx' J dfA^""^"^ {x -x',t- t') 

'Ux',t)+ilix',t) 



2n 



^^(x,t) = C'^{x,t) + — I dx' I dtT("*)(x -x',t- t') 



lheog^{x'] 



d.>A{x',t') . 



(3.6) 



Note that the integral equations incorporate the boundary conditions for the fields. The 
first set of equations implies that A{x,t) and ^u{x,t) will go to Ai„(x,t) and ^^*")(x, t), 
respectively, as t ^ — oo, and the second set implies that they will go to Aout{x,t) and 
^°"*^(x,t), respectively, as t ^ oo. 

What we shall do is solve the first set of equations for A(a;, t) and C,u{x, t) in terms of the 
in fields, and then insert this solution into the second set to find the out fields in terms of 
the in fields. 



B. Fourier decomposition 

We begin solving Eqs. (|3.5| ) by taking the time Fourier transform of both sides. Defining 

1 



A(x, uj) 



27r 
1 



dte"^'A{x,t) 
dte'^^'i^ix^t) 



(3.7) 
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and similarly for the in and out-fields, we find that 



A(x,cu) = A,„(x,a;) -^J dx' d'^^'^'^'^'^d,, ^ - 

'Ux\uj)+il{x\-uj)\ , 



l heog^{x') 



and 



Ux.u:)=i''r\x,u) + 



1 heogu{x) 1 



d^A{x,uj) . 



eoh V ~ ie — to 

In deriving these equations we made use of the fact that 



2luc 



We can derive an equation for only the field A{x, ko) by substituting from Eq. 
We find that 



(3.8) 
(3.9) 

(3.10) 
into Eq. 



A{x, uj) 



'il^-\x',u)+il^-^^ix',-u) 



2Q„ 



+ 



1 heog,{x') 



eoh V 2n^ ni-{uj + ie) 



-d.^>k{x\uj) 



(3.11) 



Our next step is to turn this into a differential equation, but before doing so we shall 
make a simplifying assumption. The field .^^*"^(x,t) is a free field which oscillates at the 
frequency Vl^, and this implies that ^\^^\x^uj) is nonzero only when uj = Vty . We are mainly 
interested in cases where the incoming light is not resonant with the medium, so we shall 
initially assume that uj ^ VL^ for u = 1^ . . . N . This implies that we can drop ^^*")(a;, uj) and 
^(«ra)t^^^ — u;) from the above equation and set e = 0. We return to the resonant case later. 

Next, we then apply the differential operator c^9^ + u;^ to both sides. This annihilates the 
in-field term and converts the integral equation into a homogeneous differential equation. 
The result is 



n?{x, Lj) 



d^A{x, uj) + uj'^A{x, uj) = . 



(3.12) 



Here, n{x, uj), the space and frequency dependent index of refraction of the medium, is given 
by 



n(x, uj) 



E 



-1/2 



(3.13) 



In this form, the equations have a rather classical appearance, and the matter operators 
no longer appear in the formulation, which gives rise to a substantial simplification. 



8 



IV. DIELECTRIC LAYER 



We now want to specialize our equations to the case of a dielectric layer with a uniform 
density of oscillators. This corresponds to the important case of a beam-splitter or mirror, 
although we make no restrictions as to the size of the layer. The medium extends from 
X = —L to X = L. Inside the medium, n{x,uj) has a value of no{uj), and outside the 
medium it has a value of 1. The solutions to Eq. ( |3.12| ) should be continuous and n~'^{x,uj) 
times their derivative should be continuous. These correspond to a continuous magnetic and 
electric field, respectively. 



A. Classical Case 



In order to find the solution of the operator equation, Eq. ( p.l2|) , we first find solutions 
of the corresponding c-number equation, which we shall denote as u{x,uj). We begin by 
dividing the line into three regions, region I for x < —L, region II for —L < x < L, and 
region III for x > L. In regions I and III, u{x,uj) satisfies 



and in region II 



c + uj u{x, uj) = 



c + uj no(c<j) u{x, uj) = 



(4.1) 



(4.2) 



where hq is the value of n{x, u) in region II. A solution incident from the left, ui{x, u), which 
satisfies the equation and has the proper continuity properties, is given by 

r e*^(^)^' + i?(tu)e-*^(^)^ in region I 
M;(x,c<;) = J 5W(cj)e^'^('^)^' + 5f^(cj)e-^'^('^)^ in region II (4.3) 
[ T{uj)e'^^^'^^ in region III , 

where k{ijj) = uj/c , and = nok{uj) . It is to be remembered that waves proportional 



to e*'^^ or e**^^ 



are propagating to the right, and those proportional to e ^'^^ or e are 



propagating to the left. 

Suppressing the frequency arguments for clarity, the coefficients in the above equation 
are given by 



R 



i{nl - 1) sin(2/tL)^„2ifcL 



T 



2n, 



^-2ikL 



(4.4) 



D D 

_ r^o(r^o + 1) (0 _ ^o(^o-l) j(.-fc)L 

- ^ e - - e 

where D = 2?t,o cos(2kL) — i{nQ + l)sin(2/tL). Note that R and T are, respectively, the 
reflection and transmission coefficients for the medium, and that + |Tp = 1. For a 
solution incident from the right we have 

T(ci;)e~*'^'^'^^^ in region I 

Ur{x,u) = <, Bi:'\uj)e'^^^'^'' + Bl''\u)e-'''^^^'' in region II (4.5) 



where k , k, R and T are as before, while Bjr'> = Bf\ and B\^^ = B^. 
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B. Asymptotic fields 



Both Ur and ui are solutions of the differential equation, Eqs. ( [4.1| ) and 
implies that they are also solutions to the corresponding integral equation 

A(x, uj] 



and this 



ic 



luj J-L 



1 



1 



d^'A{x',uj) 



(4.6) 



for particular choices of the field Ain{x,uj). We can find Ain{x,uj) for both solutions simply 
by substituting them into Eq. ( [4.6| ). We must be careful, however, because the expression 
inside the square brackets is not continuous at a; = ±L, and it is being differentiated, so 
that the discontinuities will lead to finite contributions after being integrated. One way to 
find these contributions is to consider a refractive index which is continuous, but which goes 
to the desired one as a limit. 

For example, let us suppose that n{x, u) is 1 for x < —L — S and x > L + S, is equal to 
no for —L < x < L, goes continuously from 1 to no as x goes from —L — 6 to —L, and goes 
continuously from no to 1 as x goes from L to L + S. We can then take the limit 5 — >■ 0. Let 
us examine what happens in the interval between —L — 6 and — L; the interval between L 
and L + S is similar. As 5 — we have that 



L-<5 



^iLu\x+L\/c 



-L 



dx'dr' 



L~5 



^iui\x+L\/c 



1 - 



1 ' 

v?{x, lS) 



d^>k[x\ iJ) 
d^>k[x' , w) 



d^k{x,uj] 



(4.7) 



x=-L+ 



where x = — L+ denotes the limit as x — —L from the positive direction {x = L~ is 
defined in an analogous fashion). The other limit of the integral contributes zero, due to the 
refractive index term approaching unity. Explicitly putting in the terms resulting from the 
boundaries of the medium gives 



A(x, uj) = AUx, r rfxV'^l"-"'!/^ ( 1 - 

2uj J~L+ \ nQ[ujj 



dl,A{x',cu) 



IC 



^iuj\x+L\/c I _ 



A{x,u)\ 



x=-L+ 



^iu)\x—L\l c I -j^ 



nolco") 



d.^ A{x,u))\ 



x=—L 



(4.8) 



If we now substitute ui{x,u) into this equation instead of A{x,lj), we find that 

A,„(x,^) = e^-^/% (4.9) 



and if we substitute Ur{x,uj), we find 



A,n{x,u;) = e-''^^/' . 



(4.10) 
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C. Quantum Case 



In the quantum case we have the usual expansion of a free field in terms of annihilation 
and creation operators. This leads, in the present case, to: 



dk\ 



I hceo 



^(in) ^i{kx-\k\ct) _^ ^^M^tg-i{A.'a::-|A:|ct)' 



(4.11) 



which implies that for u; > 



kin{x,Uj) 



hen 



2cAk{uj) 



(in) iojx/c _|_ " (m) -iujx/c 



(4.12) 



The results of the previous paragraph allow us to see that if A{x,uj) is given by 
A(a;, u;) = 



hen 



2cAk{iu) 



'•(in) I \ ^yiri) \ / n 



{in) 



(4.13) 



then it is a solution of Eq. (|4.8| ) with Aj„ given by Eq. ( [4.12|) . This gives us the interpolating 
field in terms of the in field. 

Our remaining task is to use the expression for the interpolating field to find the out-field 
in terms of the in-field. This can be done by substituting the expression for A(x, cj) given in 
the previous paragraph into the equation which relates the interpolating field to the out-field 

A(x, u) = Kut{x, + ^ r rfx'e-^-l^-^'l/^ f 1 - dlk{x\ uj) 

2uj J-L+ \ riQyjj) J 

k{x,uj) 



ic 
^2^ 



-iu\x+L\/ c 



1 - 



nl{uj) 



L+ 



-iLij\x—L\l c I -j^ 



x=—L~ 



which follows from Eqs. ( p.6| ). The derivation is almost identical to that of Eq. 
do not give it explicitly. Making this substitution we find that, for > 



(4.14) 



so we 



heo 



2cAk{uj) 



(r(cu)e'^(")" + i?(cj)e-*'^^"^")a^°") 



ik(u))x\p^{out) 



+ (i?(u;)e^'=(^)^ + T(u;)e )a 



ik{LLi)x\ ;^iout) 
k(ui) 



(4.15) 



The out-field can also be expressed in terms of out creation and annihilation operators, 
Aout{x,t) = 



dk\ 



hcen 



[out) ^i(kx-\k\ct) _|_ /^{oM^)^|tp-^(fca;-|fc|c^) 



ai 'e 



Taking the Fourier transform of this equation with respect to time, for > 0, gives 



(4.16) 



hen 



2cAk{uj) 



;.{0Ut) ik{uj)x I -(Ottt) ~ik(uJ)x■^ 



(4.17) 
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Comparing Eqs. ( |4.15|) and ( [4.17| ) we see that for u > 



474 - + ^H^-^H (4-18) 

«5;i) = ^MaS+^H«%) ■ (4-19) 

These equations are the solution to the scattering problem. We note that the transmis- 
sion and reflection coefficients satisfy the usual relation of |T(c(j)p+ = 1. This holds 
even in the bandgap regions, where transmission occurs via an evanescent field. Thus, the 
scattering problem is explicitly unitary, as we would expect. It is important to notice that 
unitarity holds even inside the band-gaps of the problem, indicating that the absorption 
bands simply modify the reflection and transmission coefficients, without removing photons. 
Another way to think of this, is that even when a photon is removed through virtual excita- 
tion of an atomic resonance, the photon will eventually be re-radiated in either the forward 
of backward directions. 



D. Resonances 



Let us now consider what happens when u = ±Q,^. This case was excluded from our 
earlier treatment, since it must involve some matter-operator contribution, which we have 
neglected so far. Of course, as the resonances are discrete, the frequencies involved are 
essentially a set of measure zero, lying on the upper edge of each band-gap. 

In this case, terms proportional to either ^^*"^(a;, to) or ^^*"^^(a;, —a;) will be present in Eq. 



This, in turn, means that the solution for the interpolating field given in Eq.( [4.13D 
must be modified. In particular, a term proportional to the matter fields must be added. 
This is not surprising; it is usually accepted that in a dispersive medium there must be 
absorption, and this in turn generally requires coupling to a reservoir field. However, it 
would be rather surprising to find quantum noise only occurring at discrete frequencies 
corresponding to the band edge. We have already established, in particular, that absorption 
in the band at frequencies different from the resonances simply changes the transmission 
and reflection coefficients, without adding any noise source. We will now show that even 
with the matter terms included, our earlier conclusions still hold; there are no extra noise 
sources for the out-fields. 

In order to see this, we first note that 



where: = e"*^'"* 0). This implies that: 



(in) I 



(4.20) 



(4.21) 



^1/ 



(m)t. 



-UJ 



(m)t| 



X] 



Application of the differential operator c^d^ + to Eq. ( |3.11| ), this time keeping the 
matter terms, gives: 



<9x 



d^A{x,uj) +uj'^A{x,uj) = d^Fin{x,u) , 



(4.22) 



12 



where we now include an inhomogeneous term defined as: 



= E y ^^^7^ (^(^ - + S{u^ + • (4.23) 

This equation can be solved by means of a Green's function which satisfies: 

dx — zdxG{x, x') I + uj^G{x, x') = S{x - x') , (4.24) 
\n'^{x,uj) J 

together with the proper boundary conditions. For the boundary conditions, we shall choose 
G{x, x') to have only outgoing waves at x = ±oo. Any fields produced by these matter terms 
are generated in a finite region and propagate outward, so that these boundary conditions 
are the appropriate ones. The Green's function can be found by standard techniques, and 
is given by: 

G{x,x') = 2iujcT{uj) ^^^^' uj)urix', u) [x > x'] (4.25) 
Ur{x,uj)ue{x' ,uj) [x < x'] . 



2iujcT{uj 



The solution to Eq. ( [4.22[ ) with only outgoing waves at a; = ±oo , which we shall call 
AJx, uj) is then: 



A,(x,cu) = dx'G{x,x')d^:Fin{x\uj) . (4.26) 

J — Tj 

Substitution of As(x, cj) into Eq.( p.ll ) - the integral equation for K{x^uj) - shows that it 
is, as expected, a solution with Kin{x,uj) = 0. A complete solution for general Aj„(a;, cj) can 
then be obtained by adding to A^ a solution of the homogeneous equation with the proper 
kin, as discussed in the previous section. 

As has already been noted, Fj„(x, cj) is only nonzero if = ±fij, , for some index u. 
At these values, which correspond to the upper boundary of each band-edge, the index of 
refraction vanishes inside the medium, and: 

T{uj) = e-^'^^/" (4.27) 

c — iuL 

c — iujL 



For — L < X < L we have that: 



uAx, u) = uJx, u) = ^e-^"^/" . (4.28) 

c — luL 

Therefore, for u = ±f2j, , it follows that the Green's function is constant, and As{x,uj) 
is proportional to: 

dx'd^,Fin{x\ uj) = FiniL, u) - Fini-L, uj) . (4.29) 
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Here, however, we must be careful since the source function Fin{x,uj) is proportional to 



gy{x), which is discontinuous at x = ±L , so that the right-hand side of the above equation 
is not well-defined. 

In order to resolve this difficulty, we must use the same technique as before, and consider 
a continuously changing refractive index over a small boundary region. That is, we suppose 
that n(x, uj) is 1 for x < —L — 5 and x > L + 5, is equal to for —L < x < L, goes 
continuously from 1 to as x goes from —L — 5 to —L, and goes continuously from to 1 



as X goes from L to L + 6. A similar behaviour is assumed ioi\J g^^x) . We can then take 
the limit 5^0. Let us examine what happens in the interval between —L — 5 and — L; the 
interval between L and L + 5 is similar. 

For X > L + 5 , the integral we wish to consider is then: 



/L-\-5 ^ ^ 

dx'uJx\uj)dx'Fin{x\uj) = / dx'uJx' ,uj)dxiFin{x\uj) (4.30) 
-L-5 J-L-5 



fL+S 

+ / dx'Ur{x' ,Uj)dx'Fin{x' ,u) . 



For X in other intervals, the situation is similar. Note that Ur{x',uj) is now a solution 
of the homogeneous version of Eq. ( [4.24| ), including the modified index of refraction. From 
this equation it is relatively straightforward to show that, for 6 small, 

dx'ur{x', Lu)dx'Fin{x', Uj) ^ Ur{L, u) \ Fin{L + 5,u) - Fin{L, Uj)j , (4.31) 

and this becomes an equality as 5 — . A similar relationship holds for the integral 
from —L — 5 to —L : 

/ dx'ur{x', uj)dx>Fin{x\ Uj) Ur{-L, u) (Fin{-L -5,uj)- Fin{-L, u)) . (4.32) 

Adding these contributions up, and noting that if n(x, u) = for —L < x < L, then 
Ur{—L,uj) = Ur{L,uj), wc find that: 

fL+S 

lim / dx'ur(x' ,Lj)dx'Fin(x' ,Lj) = . (4.33) 

S-^O J-L-5 

This implies that the matter operators do not contribute to the interpolating field solu- 
tion, even at the resonance frequencies which we did not consider in detail previously. Thus, 
the previous relation between the in and out operators still holds at the resonances where 
LO = . At first sight, this seems difficult to understand, since in general one would need 
to include noise operators to conserve commutation relations, and hence unitarity. However, 
this is consistent because, as can be seen from Eq.( |4.28| ), we have |T(co')p + 1/2(^7)1^ = 1, even 
when UJ = for u = 1, ..N . In summary, we reach the somewhat surprising conclusion 
that no additional noise operators are needed in the asymptotic properties of the present 
model - even at resonance. 
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V. PHOTODETECTION EXAMPLE 

Given the state of the in-field, these equations allow us to calculate the properties of 
the out-field, and hence calculate observable scattering properties. In order to see how 
this works let us consider an example. We shall find the probability that a photo-detector 
located at x, where a; > and is far from the medium, will fire at time t. At the long times 
required for propagation to this location, the fields will asymptotically become out-fields. 
The photo-detection probability is therefore proportional to 

{tn\Di-}{x,t)Dil}{x,t)\tn) = {tn\{d^Aizl{x,t)){d^Aitl{x,t))\tn) , (5.1) 

where \in) is the in state, 

AW(x,t) = / dk^^^^at\^^'^-\'\-^ , (5.2) 

and A.oJ^{x,t) = (Ao^](x, t))''". Now let f{k) be a function which is zero if < . The 
Fourier transform of f{k) is closely related to the shape of the pulse which is being sent into 
the medium. Define 

alif] = I dkfik)iat^y , (5.3) 

and let 

I in) = exp(ai„[/] - aj„[/])|0)i„ . (5.4) 

This is a coherent state composed of wave packets with the intensity of the field and the 
shape of the wave packet determined by f{k). For this state the correlation function in Eq. 
( ^.1| ) is given by 

(5.5) 

where we have used Eq. ( [4.19| ) to relate the in and out operators, and we have explicitly 
indicated the k dependence of the transmission coefficient. 

As expected, this equation demonstrates explicitly that photodetection rates are sup- 
pressed for frequency components that correspond to the dielectric absorption bands, where 
T{uj) — > 0. At these frequencies, the predominant effect will be a strong reflection, with no 
photodetection occurring at the detector location on the other side of the mirror. 

VI. CONCLUSION 

We have presented an analysis of a quantized electromagnetic wave scattering off a linear, 
dispersive medium of finite extent. The medium consists of harmonic oscillators whose 
energy-level spacings can be chosen to match those of an actual medium in the spirit of the 
classical Sellmeir expansion. What emerges is a relation between the in and out fields which 
is most simply stated in terms of their annihilation operators. The overall results are exact, 
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and simply expressed in terms of linear transmission and reflection coefficients. The medium 
has both transmission and absorption bands. Since the model is a full quantum-field version 
of the widely used Lorenz model that leads to the Sellmeir expansion, it has a wide area of 
applicability to realistic dielectric media with a variety of dispersion relations. 

The final results are very similar to the classical expressions which relate the amplitudes 
of the incoming and outgoing waves. The transmission and reflection coefficients which one 
finds from the quantum and classical analyses are identical, as they should be, given that 
the model is a linear one, and must reduce to the classical theory in the correspondence 
limit. Expressions such as those appearing in Eq. ( 4.191 ) are often used in quantum optics. 



and we believe it is useful to see how they emerge from the underlying scattering theory. 

An important feature of the theory treated here is that it has absorption bands, without 
any corresponding noise terms in the field equation. This is due to the dielectric model 
used here, in which the dielectric constant is always either purely real or purely imaginary. 
In this type of model, all photons absorbed are re-emitted. Thus, absorption simply re- 
sults in a strong reflection, with purely evanescent fields inside the dielectric. In a related 
phenomenological treatment , one finds similar behaviour: if the dielectric constant is 
always either purely real or purely imaginary, it is possible to have dispersion without any 
additional noise terms. 

The reason for this is that the scattering terms alone are sufficient to ensure that the 
input-output relations remain unitary, with no change in the commutators. Thus, in the 
present model, there is no need for any additional source terms. The theory therefore 
provides a justification for the use of simple input-output relations to describe idealized 
dielectric or metallic mirrors, even when the dielectric response is dispersive. However, for 
realistic media it is generally the case that absorption can also occur even in the transmission 
bands. Treating this would require the use of more sophisticated models, including a complex 
refractive index. 
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